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TECHNICAL NOTE 2620 

PRINCIPLE AND APPLICATION OF COMPLEMENTARY ENERGY METHOD 
FOR THIN HOMOGENEOUS AND SANDWICH PLATES AND 
SHELLS WITH FINITE DEFLECTIONS 
By Chi -Teh Wang 

SUMMARY 


The principle of complementary energy in the nonlinear elasticity 
theory is shown to he derivable from the principle of potential energy 
hy a Legendre type of transfonnation. In particular, the expression of 
the complementary energy is derived for homogeneous and sandwich plates 
and shells with large deflections. By the method of complementary 
energy, the stress-strain relations are derived for homogeneous shells, 
sandwich plates, and sandwich shells. Without the use of this method 
much lengthier calculations would he necessary. 


INTRODUCTION 


In the theory of elasticity, the most important variational 
principle is perhaps the principle of potential energy, which states 
that of all displacements satisfying given houndary conditions those 
that satisfy the equillhrium conditions make the potential energy a 
stationary value. For stable equilibrium, the stationary value may be 
shown to be a mi ni Tmim (reference 1). The potential energy is defined 
as the difference between the strain energy eind the potential or virtual 
work which the surface stresses do over that portion of the surface on 
which the surface stresses eire prescribed. This principle is capable of 
general application as it holds true no matter what the law connecting 
load and deformation may be (reference 2) . With the relationships 
between stresses, strains, and displacements known, the differential 
equations defining the equilibrium conditions may be derived from the 
variational principle by the methods of the calculus of variations.- 

The principle of potential energy was obtained by comparing the 
strain energy U of the equilibrium state, characterized by displace- 
ments u, V, and w, wl-th the strain energy U + of a nei^boring 
displacement state u + Au, v + Av, and w + Aw. A corresponding 
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variational principle may "be derived by varying tbe stresses rather 
than the displacements. This results in the so-called principle of 
cong)lementary energy which states that of all stress states satisfying 
the conditions of equllibri\im in the interior and on that portion of 
the surface on which the surface forces are prescribed the actual state 
of stress is such that the complementary energy is a stationary value. 
AgEiin, for stable equilibrium, this stationary value is a minimum. 

For a material that obeys Hooke's law and has linear strain-displacement 
relationships, the complementary energy is the difference between the 
strain energy and the work which the surface stresses do over that 
portion of the surface where the displacements are prescribed. However, 
for nonlinear problems, it seems that the expression of the complementary 
energy has not been given previously. 

When the equilibrium equations are known in terms of stresses, the 
stress-displacement or stress-strain relations may be derived from the 
principle of conplementary energy by methods of the calculus of variation. 
In some elasticity problems this approach has been found to be convenient. 
For example, in the linear case Trefftz (reference 3 ) lias used the 
method for the derivation of the stress -displacement relationships in 
the case of thin homogeneous shells and Eeissner has recently derived 
the stress-displacement relationships of sandwich plates (reference 4) 
and shells (reference 5 ) hy the same method. Without resource to this 
method these derivations may require much lengthier calciilations . 

In this report the principle of complementary energy is derived 
for thin plates and shells with large deflections. The plates and shells 
may be either homogeneous or of sandwich-type construction. The appli- 
cation of the principle to the derivation of the stress-displacement 
relations for homogeneous and sandwich shells as well as sandwich plates 
is given. 

This work was conducted at the Daniel Guggenheim School of Aeronautics, 
College of Engineering, New York University, under the sponsorship and 
with the financial assistance of the National Advisory Committee for 
Aeronautics. The author is Indebted to Professor K. 0 . Friedrichs for 
his discussion of the problem and to Dr. G. V. K. Rao for his helpful 
assistance . 


SYMBOLS 


A total area of domain 


C 


boundary of domain 
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D flexural rigidity of plates and shells (Eh^/l 2 (l - V^)) 

e effective transverse normal strain for sandwich plates and 

shells 

E Young’s modulus of elasticity 

G modulus of rigidity in shear 

h thickness of homogeneous plates and shells; also thickness 

of core layer in the case of sandwich-type construction 

I potential energy 

J complementary energy 

Mxj%-,Mxy resultant bending and twisting moments in plates and shells; 
unit in moment per tinit length 

Mnv.MTiy X- and y-conrponents of the resultant moment on boundary 

Wx,Ny,Wxy stress resultant in middle plane of plate or in middle surface 
of shell; tmlt in force per unit length 

p external lateral pressure 

q difference in p acting on upper face and lower face of 

sandwich plates and shells ^(Pu - 

shearing stress resultant in plate or shelly normal to middle 
surface, unit in force per unit length 

Qn resultant Q on boundary 

R radius of curvature at any point in shell 

t thickness of face layer of sandwich shell 

U work done by stress resultants due to large deflection 

V strain energy 

W virtual work done by external forces and moments 

x,y,z coordinates; y is cvcrvilinear coordinate in direction of 

circumference in the case of a shell 
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Xn,Yn 

u,v,w 

Px^Py 

7 

6 

V 

V 

Xl, ^l|-i ^ 

cfz 

**XJ "V 
X 

Sutscripts 

c 

d 

f 

I 

s 


o 

u 

in 


X- and y-components of stress resultant acting 
on ‘boundary 

displacements in x-, y-, and z-directions, 
respectively 

components of change of shape of normal to middle 
surface of plate or shell 

shear strain 

first variation 

strain in x- and y-directions, respectively 
Poisson's ratio 
Lagranglan multipliers 

component of transverse normal stress in core 
layer 

transverse shearing stress in core layer 
arbitrary parameters 

change of curvature at any point in shell 

core layer of sandwich plate or shell 

portion of boundary over which displacements are 
prescribed 

face layer of sandwich plate or shell 

lower face of sandwich plate or shell 

portion of boundary over which stresses are 
prescribed 

given values on boundary 

upper face of sandwich plate or shell 

values of stress congponents at middle surface 
z = 0 
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E3CPRESSI0N OF COMPLEMEHTAEY EMERGY FOR 
THIN HOMOCMEOUS PLATES WITH LARGE DEFLECTIONS 


The principle of coraplementary energy can he derived from the 
principle of potential energy 1^ means of a Legendre type of transfor- 
mation. Since the principle of potential energy is also valid for 
nonlinear strain-displacement relationships, the potential energy 

I = V - W 


is stationary or 51 = 0 when I is varied with respect to the 
displacement components u, v, and w. In the above expression V 
is the strain energy and W is the virtual work. In the case of a 
thin elastic plate with large deflection the strain energy consists of 
two paxts, one due to the stretching of the plate and the other due to 
bending. Consequently, the potential energy may be written as follows 




Nx^ + Ny^ - 2VNx% + 2(1 + V)Nxy"' dx dy + 


+ V)Nxy^ 


2(1 


- 2VMx% + 2(1 + V)Mxy" 


dx dy - 


1 


pw dx dy - 

A JCa 


XnU + YnV + fXn ^ + Yn ^ + 


M „ ^W 


ds 


( 1 ) 


where A is the area of the plate and Cs is that portion of the 

boundary where the surface forces are prescribed; I has a stationary 
value when it is varied with respect to u, v, and w, subject to the 
following stress -displacement relations: 
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Wtt = 


Wxy - 


Eh 

r 

Bn 

+ 

fBw\^ 

4“ V 

By 

1 - v2 

Bx 

2 ' 


T V 

By 

Eh 

Bv 



+ V 

Bu 

1 - 

By 


;Byj 

Bx 

Eh 

1 

fBz H 

1 Bu . Bw 

Bw^ 

2(1 + 

V)' 

bx 

By 

Bx 

ByJ 


Mx = + V 


= + V 

W^ 3x2/ 

Sx ^ 


“y 


Mxy = D(1 - v) 


= D + ^'1 

^\Sx^ Sy^/ 


> 


(2) 


^ = D |-(& + ^ 

^y^/ 


Expressions ( 2 ) are the usual stress-displacement relations when the 
deflection of the plate is large compared, with its thickness but is 
still smal l enough to allow the use of simplified formulas for curvature. 
The boundary conditions are: On that portion of the surface where the 

surface forces are not prescribed, boundary displacements u, t, and 
w are equal to the given values Uq, Vq, and Wq. 

The principle of complementary energy may be obtained directly by 
applying a Legendre type, transformation to the principle of potential 
energy. This procedure is known as "Friedrich’s Method" (reference 6) 
in the calcnilus of variation. Let ^ 

be the Lagrangian multipliers and px.> similar 
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parameters. Then the variation jirohlem., equation (l), is equivalent 
to varying the integral H with respect to Hxy^ 

u, V, and w independently, where 



- 2vHxKy + 2(1 + v)Kxy‘ 


G 


dx dy + 


2(1 


Mx^ + My2 - 2VMxMy + 2(l + V)Mxy2 dx dy 


+ V — + 

hy 2\^ 


^ f - ^ ^ Kv) ' ^ ^ i ^ i(i)!l 


^ dx dy + 


L dx dy + 




N- 


xy 


Eh ^ ^ ^ ^ 

2(1 + v) ^ ^ hjj 


dx dy + 


k02! 


I 


^w . 


P4x - D + V 




+ V 


dx dy + 


dx dy + 


Mjjy - D(1 - V) 


■ dx dy - FP pw dx dy - 

ox ^ vA/a 


jT jx^u + Y^v + ^ ^ + Qn^w + Mxn g + Myn ^ 

^ j^l(u - Uo) + - '^o) + f^(w - Wq) + 


ds + 
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The condition 5H = 0 leads to: 
6E = 0 


-i 


mjr - 


Eh 

pg v)Ng ^ + J 
L Eh 3 Lg-v2)D 3 

[%■ - r_^ . Xglaw 

pi - v2)d ^ ^ pi - v)D J ^ 

1 _ ^ 2(1 + v)J 


VN-y \ /^V “■ ^ 

- 

L(i- 


6Ny + 


+ A-A 


s ^ s Eh 


1 


. (X2 + VXi) &T + 

8x 2(1 + v) 8y 1 - _ 


.r.^r ■ ^ , 

1 - v2 (1 + v) 8x 8y 


Eh(l2 + VX-i) ^ ^ 


-V^ ^y2 8x^ 1_V- 


^ Eh(Xj^ + VX2) ^ EhX^ 


8y 2(1 + v) 


8w 


EhX3 

3 Eh(A^ + vXi) 

5y 

8x 

2(1 + v) 

8y 1 - v2 


- p ^ 5wj dx dy - 


I 


1 -^(Xi + VX2)SU 1 -^5^ 

8xl-v^ 8y2(l + v) 


A . A _S_(X2 + VXi)hv 

8x 2(1 + v) 8y 1 _ v^ 


A ^h (Xi + VX2)^ 8w + A _s^3_ Ja 6w 

8x 2 _ v^ ^ ^ 2(1 + V ) 8x 
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^ ^ S Eh 

Sx 2(1 + v) Sy ^ 1 - 


(A.1 + VA.2) 


^ 5w 


+ 


D(X4 + vX5)6 ^ + a d . ( . 1 - - - ^) 
^ Sy 2 








Iw 

•y 


•^8 ^ D(X5 + VX4)5 ^j + 


|L 6w D ^ Sw D ^ 


Sx 


Bx 


^^x^ ay^ 




^\ax^ 


^dx dy - 



Xn5u. + Yj^Bv + 



^ + Yn ^ + lin^Sw 




ds + 


J^ I-I28V + iJ^Sw + ii.45 ^ + M^S ds 


w 


where 


„ a2 ^2 
ax^ ay^ 

In carrying out the above C60.ciilation, integration hy parts has 
been used, an example of which is 


1 

1 


Eh 


1 - V 

Eh 

5x 1 _ -V 


* VX2)^ ^ 

^ '^'^ax 2(1 + V) Sy 


dx dy = 


-(Xl + vX2)6u + 


EhX3 


6 u 


F 


a Eh ^ ^ 

^ “ 


ay 2(1 + V) 

5u dx dy 


dx dy + 


1 - V 


ay 2(1 + V) 


(5) 
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According to the veil -known rules of the calcul\is 
Euler equations of the variational prohlem are 

of variations. 

■ "“y + XI - 0 

Eh 

(6) 

Nv - VEx 
^ + \2 - 0 
Eh 

(7) 

2(1 + v)Exy 
+ X 3 = 0 

Eh 

(8) 

0 

II 

r< 

1 1 

(9) 

My — 

^ + Xr - 0 

(1 - v 2 )d 

( 10 ) 

2Mvy 

X 6 = 0 

(1 - V)D 

( 11 ) 

Eh . ,n ^ . ^ ^3 

3x 1 . ^ * Sy 2(1 - V) - ° 

( 12 ) 

a EM 3 ^ a Eh , , 

3x 2{H- V) ay 1 + 

( 13 ) 


„ . Eh(Xi + VX2) ^2^ EM-3 , 

DV^ + — pr + ^ 

1 - ^ (1 + V) ^ ^y 

Eh(X2 + V^l) ^ ^ Eh(Xi + VAg) ^ ^ EhX3 ~ | ^ 

1 - ^ 1 -v^ ^ 2(1 + 

^li_ ^^3 ^ i. Eh(Xi -h VX2)| _ ^ ^ ^ 

Sy Sx 2(1 + v) ^ 1 - v^ J 


( 1 ^) 
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With the Lagrangian multipliers Xp, . . . X 6 determined by equations ( 6 ) 

to ( 11 ), equations ( 12 } to (l4) become identical with the equilibrium 
equationsj namely 


^xy 


0 


(12a) 


SWxy 2®y 
hx. hy 


0 


( 13 a) 


DWItf + W- 




+ 2W- 


xy 


hx. ^ 


+ H- 




p = 0 


(llta) 


In writing down equation (l^) the relations (l2a) and (l3a) have been 
included. 


Next, substitute the relations ( 6 ) to (ll) into the second group of 
surface Integrals in equation (4). By Gauss’ theorem these surface 
integrals may be transformed into line integrals as follows: 



+ VX2)6u + 


^ EhX 36 u 
^y 2(1 + V) 


dx dy = / ^ f ^ UxSu + — W^ybuy dx dy 




6 u ds 


-I 


!_ 6v + A ^-(Xe + VXi)Sv 

^ 2(1 + v) ^ 1 - 


dx dy 


-L 


8 u ds 
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-I 


1 frl + VX2) 52 5w + A 5s 6w 


1 - 


^ 'by 2(1 + v) Bx 


dx dy = / Xji — 6w ds 
'-'C ^ 


-1 


®“*-3 52 8W+A “ 


bx. 2(1 + V) ^ 


(X2 + VXi)— Sw 
Sy 1 _ v2 ^ 


dx dy = J' Yn — 6w ds 


bfv 

'C ^ 


.,•3 


-£ 


b / \ bw S D(1 - V) -> -e bw 

^ D(X 4 + Vl 5)5 — + — ^ ^ -- ^ 6 ^ — 
bx. dx ay 2 


Sx 


dx dy = / Mhx 5 ds 
U n 


hx. 


-f£ 


A Xg5 + :^ D(X5 + VX4)8 




by by 


dx dy 


= MjiyS 


bw 


ds 


- /TIA 8w D + 5!»\ + A 6» D A/k!^ + 5!«\ ax ay - r Sn6w da 

'-'A ^ ^\^x^ ^y^/ ^ ^\^x2 byy 


Combiiiiiig these integrals with the line integrals in equation 
the condition that the line integrals vanish requires that 


^ll = -Xn 


= -Yn 


^^3 



+ 



+ 



= -Mnx 

^5 ~ “^ny ( 15 ) 


on C(i. 
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Substituting relations ( 6 ) to (I5) into equation ( 3 ), carrying out 
integration "by parts, and applying Gauss’ theorem to transform surface 
Integrals into line integrals wherever possible, one obtains finally 

H = - - 2 vHx% + 2(1 + v)Wxy^ dx dy - 


1 

2(1 - v 2 )D 



+ My2 


2-AMx% + 



dx dy - 



Mnx ~ + Mjjy — 

ox oy 


ds 


(16) 


The complementary energy J is defined as negative H or 


J=-H=V+U-W 


(IT) 


where V is the strain energy, U represents the work done by the 
forces in the middle plane of the plate due to large deflection, and W 
is the work done by the surface stresses over that portion of the 
boundary where the displacements are prescribed. 
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APPLICATION OF PRINCIPLE OF COMPLEMENTARY ENERGY 
Thin Homogeneoiis Cylindrical Shells 


Assume that the complementary energy for the circular cylindrical 
shell may he ■written in the same form as in "the case of a flat plate. 
With the equilibrium equations derived., the stress-displacement relations 
may he obtained from the principle of - con^plementary energy as follows; 

Consider the case where w is large in comparison "with u and "v. 
Retaining the terms consisting of products of derivatives of w times 
the stress res'ultants and co'uples, the equilihri'um equations for an 
element of a cylindrical shell with arhitrary cross section (fig. l) 
can he put in the following form: 


^Ny ^Nyv 

— + ±i. + =0 


(18) 


SNy ^ SOljjy Njjy ^ Q 

Sy Sx R Sx R 


(19) 




Sx ^y ^xV"^^ 8x dy, 








^ Bw Ny 

^ + P = 0 




8Mx ^ 8Mjjy ^ My ^ 

hx. dy R ^ 


— - 0 ^ = 0 


(20) 

(21) 


2 m- 


Bx ^y R dx 


( 22 ) 
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wliere x is in the direction of the generator and y is in the direction 
of the circumference. These equations may he obtained following a 
similar derivation as given hy Timoshenko (reference 7 ) with slight 
modification. Let X3, llj., and X5 he Lagrangian multipliers. 

5y varying the stress resultants and couples independently in the 
con5)lementary energy equation subject to the equations of equilibrium, 
one finds that the process is equivalent to formulating a variational 
integral L and setting its first variation equal to zeroj namely. 




|Wx^ + - 2vWx% + 2(1 + v)Nxy^J dx dy + 


"grWlE-- + - 2-\M3^4y + 2(1 + dx dy 


€ 




+ 2Nxy ^ 




hw hw 









(I 


dx dy + 






R ^ Ry 




1 4 . J. ^ L ^ 

8y 

A(h„ |£ t n„ St |£ -I. Sit 4- pi 

3 yy ^ ^ ^ dyy R ^ R _j 


^SMx % ^w \ ^ 




f F 

^CdL 


XnU + Y^v + (Xn ^ + Yy ^ + Qn]w + MjjxPx + Mny^y | ds = 0 






(23) 


where the surface integrals are to be extended over the entire cylindriceil 
surface, 3x angle through which Mjjx turns, and Py is the 

angle through which Mj^y turns. 
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Carrying out the first variation, regrouping the terms, and 
eliminating the variations of derivatives hy integration hy parts, one 
obtains : 


6 L 


- /T - v% _ ^ ^ ^ 

Eh ~ Sx! 


8N-y + 


Ey - VRjj. ^ ^ 


Eh 


^ R ' 2W 3y 


^^3 , ^ 

^ R ^ 


6Ny + 


\ 

2(1 + y) ^ ^X- 2 ^ X2 ^ ^ ^ ^ ^ ^ 

i Eh ^ ~ ^ ^ R^x^Sy^^ 2iy ^x_ 


-X4 - 




8Wxy + 


Mx - "W^y 

^Xlj. 

6Mx + 

- \Mx 

_ ^ ^ 

(1 - v2)d 

Sx 

|ji - v2)d 

^ R ^ 

2Mxy 

^X4 


^5 ^ ow 
— — 8Mxy 

R 

► dx dy + 

(1 - V)D 

Sy 

3x 


5My 


rr 7 bXi5Nx ^ SXl5Rxy\ ^ /^X25Rxy ^ ^ 

/SXgBQx SX 28 ^\ /SXij.5Mx ^X,l(.S4xy\ 

\Sx Sy/\^ 

/BX5^ ^ BX5^\ ^ ^ _ r Csx^ + + 

\ Sx hx /_ '-'C^ *- 

8Xn + — SYn + eOnV + Px^nx + Py^vJ ° 
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Reducing the second surface integral in equation { 2 k) to a line integral 
hy Gauss' theorem, the final hoixndary integral of equation { 2 h) becomes 


J + A25Yn + ^ 8Yn + &Qn^ + 


ds - 


/F 

^CdL 


uSXn ■+ v5Yn + wfe ^ + 5qJ + Px^M^x + PySM^y 


ds (25) 


The first line integral vanishes on the part of the boundary where surface 
stresses are prescribed. It therefore has a nonzero value only on that 
portion of the boundary where the displacements are prescribed. 

Since expression ( 25 ) must be zero, one obtains on the boundary 


Xj[ = u 

X2 = V 

X3 = w 

^4 = Px 

45 = Py 


(26) 


As equation ( 24 ) also holds for any part of the cylinder, if the 
boundary displacements referring to this part are identified with the 
displacement occurring in the actual solution of the problem, it follows 
that the Lagrangian multipliers throughout the shell are related to the 
generalized displacements in the interior of the boundary. It follows 
that equations (26) also hold true now in the Interior. 

Usirg equations (26), the Euler equations of 5 L = 0 give the 
following stress -displacement relations of the cylindrical shells: 
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Kx - VMy _ ^ ^ l/8w\^ 
Ell 2 \8x/ 




Hy - VNx _ ^ V ^ l/^\2 u ^ 

8y R 2\8y/ R 8x 

2(1 + v)Rxy _ ^ ^ ^ ^ ^ ^ V 8w 

Eh ^ By Bx By R Bx 


_ T ^ ^ ^ 

^ R By R 8x 

Mx - v% BPx 

(1 - v2)D 8x 




(27) 


- VMx ^Py Px 

(1 - v 2)D ^ ^ 

^xy _ ^Px ^ ^ Py Bw 

(1 - v)d 8y 8x R 8x 


By means of Hooke’s Law and the relations between moments and 
changes of curvature, equations (27) can he reduced to give 
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ex = (Nx - V%)/Eli 
- ^ 

2ySx/ 

6y = (% - VNx)/Eh 

= - — + l/^\ ^ u ^ 

^ R 2\^/ R Bx 

7xy = [2(1 +V)Nj/Eh 

^x hy hx. hy R Sx 
= _(Mx - vMy)/(l - v2)D 

" 5x2 

Xy = .(My - VMx)/(l - v2)D 

_ 1 ^ ^ 52y 5 /w 5w\ ^5 w\2 

R 5y 5y2 5y\R 5x/ 'K\^xj 

Xxy = -Mxy/(1 - v)D 

= -i- ^ _ i ^ 1 5w 5v 

2R 5x 5x 5y 2 5x\R 5x/ 2R 5x 5y 


It is interesting to note that the above expressions contain terms 
which have not been taken into account previously. Nor is it obvious 
from a study of these expressions that these terms are necessarily s m a ll 
compared, with those retained. 
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In the above calculation^ the expression for the complementary 
energy ban been assumed to he of the same form as in the case of a thin 
flat plate. This fact will now he verified in the following manner. 

The principle of potential energy is known to he applicable in this 
case. If incorrect stress-displacanent relations are used the correct 
equilihrium equations will not he obtained by means of the principle of 
potential energy. Suppose the Incorrect expression of the complementary 
energy was usedj the stress-strain relations obtained natxirally would 
not be correct. With these incorrect, relations, the first variation 
of the potentieil energy will not lead to the conrect equilibrliim equations. 
However, in the present case, when equations (28) are used, the first 
variation of the potential energy leads directly to the equilibrium 
equations ( 18 ) to ( 22 ), thus confirming the fact that the expression of 
conplementary energy used is correct. 


Cylindrical Shells and Flat Plates of Sandwich Construction 

As in the case of homogeneous thin shells, the stress-displacement 
relations for sandwich shells may also be obtained from the principle 
of complementary energy. A sandwich shell consists of a core layer 
surrounded by two face layers. The face layers are treated like thin 
shells of thickness t having negligible bending stiffness about their 
own middle surface. The loads applied to these face shells are of two 
kinds: Exrbernal loads and loads caused by the stresses in the core 
layer. The core layer of thickness h is assumed to behave like a 
three-dimensional elastic contlnuvtm in which those stresses which are 
parallel to the faces are negligible compared with the transverse shear 
and normal stresses. (See fig. 2.) 

Assuming a large deflection, the equilibrium differential equations 
for the upper face layers can easily be shown to be 


^^xu 

8x 


1 + 


h + t 


2R 


^%xu ^ Nyu 8wu 
8y R 



h + ti 


2R 


= 0 


(29) 


^IJxyu / ^ h + t\ 
8x I 2R / 



= 0 


(30) 
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Nyu S 
H Sx 


Nxu|l + 


h + t\ 


&W| 




2R/Sx 






N- 




xyu 


u 


%U Swji 


Sx ^ h + t ^ 
2R 


L ^ + t. 

Pu 1 + — r^l + ll + 


2R 




"yui 


h + t \ . _ . h 


2R 


+ T. 


XU 




1 + 


±l\ 

2R ) 


( 31 ) 


Similarly, for 1116 lower face layer. Hie equilibrium equalious are 


2R/ Sy R^x \ 2R/ 


. iL±l\ . 521 5^1 . Tyj/i - = 0 (33) 

Sx\ 2R J by R^ \ 2 R/ 


R 


S C ( h + t\ „ Sw7~l 




ST'" 


”yi 






PZ 




SW7 

-.z 5ji|i - 


bw-j 

2R / ^ ^ ° 


( 3 ^^) 


where Ihe subscriplB u and I denote the upper emd lower surfaces, 
respectively . 
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Under the assumption of negligible face -parallel core stresses, 
the equilibrixDu equations for the core layer are (reference 5) 



Let the values of the three stress components at the middle surface 

(z = 0) be designated by the subscript m. Integration of equations ( 35 ) 

"to (37) results in 



Let Ox and Qy be the resultant transverse shear forces. Then 




h+t 

'/ f ' 5) = 

2 


h+t 


r 2 (h + t) ”*yni 

Ty dZ = ^ 

^ , fh + tV 

2 ^ 2R / 


ihl) 


{h2) 
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From equations (28), ( 39 ), (^ 1 ), and ( 42 ) the following relations may 
he obtained: 


1 + 


2R y'xu - ^ 


= 0 




R 


> m 


h + t 


(1 - - ^)- 


x2 


h + t 




Ox 


- %■ 


t. 4 ~ 

For thin sandwich shells « 1, In such cases a combination 

2 R 

of equations (4o), (4l), and (42) gives 


(1 * - (. - ^ (44) 


In view of the fact that all face-parallel core stresses are 
neglected, the face-parallel stress resultants and couples of the 
composite shell are due to the stresses in the face layers only and 
may be obtained as follows : 
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% - Nyu + Ny 7 


= Wxyu(l + + Nxyl (1 - 


Nyx - Nyxu + ^yxZ 


Mx * * ^) * %2 (1 - 

My = H^(-V + %z) 

“xy = 2 ^ jj®xyu + ®xyj ^ 

%x = ^ 2 ' - ^-^yxu + %xl) 

P - Pu(l + + PZ (1 - 

““ - 5^) + Oxl (1 - ~)_ 


Since Uyxu = ^xyu^ ^yxZ = and « 1 , one may vrite 

®^xy = ^yx ^xy = %x* 



4R 
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Define the deformations of the composite shell to he w = + ■w’2)/2 

and e = (wji - 'W2)/(h + t) -vdiere w represents the effective transverse 

'deflection of the middle surface and e represents the effective trans- 
verse normal strain for the con^osite shell. From these definitions, Wu 

and -w^ may he written in terms of w and e as follows: 


w^ = w + 


(h + t)e 


Ml = V 


(h + t)e 
2 


(h6) 


With equations (^3), (4U), (^5), and (46), the following equations 
may he obtained hy carrying out addition as well as subtraction of 
equations (29) and (32) and (30) and (33)- 


SNx ^ ^Nxy ^ ^ ^ ^ ^ _ Q 
^ 6y E 6x R 6x 

^xy ^ ^^xy ^ ^ 6e %• _ Q 

^ ^ R 6x R 6x R, 

^xy Q 

Sx 5y R Sx ^ 

2Mxy Mxy 6w (h + t)%xy 

6x 6y R Bx 6x ^ 


0 


(47) 

(48) 

(49) 

(50) 


From equations (31) and (34) the following two relations may he derived: 




--I 


+ 



6w 

6x 



8x 





+ 


p + 


/ ^Qx 

\8x 8y / 


% ^ 
R 8y 




= 0 


(51) 
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15 _lD£ + h, 

4 Sxy ^ 






3^ t S„ a»\ (ht_^ % , 0 


(h + tXozm + 9 ) - ^ j, 


( 52 ) 


When the radluB of curvature R hecomes infinite, equations (47) 
to ( 52 ) reduce to the same equations obtained hy Reissner (references 8 
and 9) . I'Then the effective transverse normal strain e for the composite 
shell is neglected, equations (47) to (51) reduce to the^ usual equations 
of force and moment equilihrium- Equation (52) has no counterpaort in 
the theoiy of homogeneous shells, in the sense that the corresponding 
equation for the homogeneous shell contains Information that is not of 
practical interest and is therefore never formulated. This equation 
gives the local change of thickness of the shell caused directly hy the 
external loads hy way of the nonlinear terms having stress resultants 
and couples as factors. 


Denote the properties of face layers end the core hy the suhscripts 
f and c, respectively. The strain energy of the face layers and the 
core may he written as (reference 5) 


Vf = i + 2(1 + Vf)Nxy^] X 

- 2 v 0 l 3 ^y + 2(1 + Vf)Mxy^ I dx dy 

t(h + t)%!f L ^ 




^ h + t 

(h + t)Gc Ec 


2 1 f^Qx ^Qy) 


^zm + 


12\8x 8y / J 


>dx dy 


The work done hy the stress resultants due to large deflection is 


" = !// 




(8wu\2 


'xu 


k8x 


+ 2N 


8wu 8wu 


xyu 


8x hy 


+ R 




f^uj ^ 

fl + 

W / _ 

V 2R / 


fl - ^ 

W / _. 

\ 2R / 


dx dy + 
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With the aid of eqiiations (4^5) and (46), U hecomes 


0 = i + 2H„ Js Sm + 


2jj r\^l 


dx dy + 


(h + t )2 


8 


■If 


Wxl — 1 + 2N. 




xy 


^ ^ ^ / 8 e\ ^ 

8x 6y \^8yi 


dx dy 


F 


8 e 8 w / 8 e 8 w 8 e 8 w\ 


Mx ^ ^ + Mxvl— — + :i:xrzi + Mv — — 
8x 8x .. \8x 8y 8y 8x, ^ 


8e 8w 
8 y 8 y 


dx dy 


The work done hy the surface forces over that portion of the surface 
where the displacements are prescribed is 



where the last two integrals are the work done by the shear stress 
resultant in the core. The first one of these two is the work done 
due to the displacement of the middle surface and the second, due to 
compression of the core. 


Define the displacement of the middle surface and the resiiltant 
load on the boundary of the con^josite shell to be as follows: 


u = |(uu + Ml) 


Y = l(vu + Vj) 


3x = 


(-Uu + Uj) 


h + t 
Xn = Xnu + Xn2 


^nx - 'Z (~Xjiu + Xqj) 


Py = (-^u + V2) 

h + t 


Xn = Xnu + Xnl 
h + t 


Mny - 


(-Xnu + Xnz) 


2 


2 
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The work W "becomes 


^Cd 


Xnu + Ynv + /x^ ^ + Yn ^ + QnV + 


e + 


■^WPx ^ + Mny ^1''^ “ (^'*nx ^ + Mijy ^ 

IiL±_^L ^ + Y^ |^]e + + ^)qh 

4 ^ 8 x 837 12 Ec ^ / 

Introducing the Lagrangian multipliers to Xg and cariying 

out the first variation with respect to Nx, Ky, ^x^ My, Mxy, 

Qy-j and 02 m independently, the relation 

5L = 8 Vf + 8 Vc + 6 U + SV7 + S /T<jxif— + ^ + 

JJ\ W ^ R Sx R Sx/ 


^ 2 ^ 


*^xy Nxy ^ ^ Mxy ^ ^ 

8x ^ R^ R8 xR/ 


^ ^ + 
® 8 x \ 8 x/ 


8 


^ + i- 

8 x/ 8 x 


8 e' 


••• 


(-5)-sf'S)-s(”-S)-5(^5 





Xk 


^ ^xy^ ^ ^ (h + t)^ ®y 8 e 

8 x 8 yR 8 x 4 r 8 x 


+ X^ 


SMxy SMy 
etc 8 y 


8 w^ 

'by) 


+ M;xy — 1 


^ 3 w ^ ,(h * t)£ ^ J ^ ^[% ^ Sw 

R 8 x 4 R 8 x _j |_R ^ \ 8 x 


8 e 


Ry I - (h + t) (Otthi + 1 ) 
by/ 


■(•■fi 


+ Qj, 


by) 


R 


dx dy = 0 
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leads to tlie conditions = u, ^ t ^3 ~ A 4 - ^ Py# 

and Xg = -e and the following Euler equations 


Nx - Vf^^y _ Su , 1 


2tE 


bx 2 


4 . (h + t)^ /de^^ 

bxj 4 Idx, 


% - ^ _ w + 1 

2tEf ^ R 2 


+ (h + 

^y/ ^ V 


u ^ ^ (h + t)^ Px 5e 
R bx k R 5x 


2(1 + '^f)Rxy _ ^ ^ bv (h + t)^ ^ ^ + Z ^ (h + t)^ Py Se 

2tEf by bx bx by 4 SxSyRSx 4 rSx 


Mx • 

- ^f^y 

_ 5px 

_ ^ 

Se 

t(h + 

t) 2 Ef /2 

Sx 

5x 

5x 

% ■ 

- VfMx 

. ^Py 

bw 

ae ^ 

t(h + 

t)%f /2 

Sy 

by 

ay 

4(1 + 

Vf)Mjy _ 

PPx ^ 

apy 

aw 

t(h + 

t) 2 Ef 

dy 

ax 

ax 


_ ik 




^ _ Zi ^ + e 
bx R 3x R 


(h + t)Gc dx\3x 

^ _ (h + t) bl^ 

(h + t)G(, 12Ej, 


bw be 
' by ~ by 


bw 


w dw 


iSv + e- S”- = ^ + i + 5 §* 


It-pL s - e . 

4 R 


> (53) 
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In the case of a flat sandwich plate (R — ) <») the stress-strain 
relations (53) reduce to exactly those obtained hy Eeissner (reference 8). 
In the case of homogeneous shells (Gc = Eq — ^ e — ^ O) equations (53) 

then reduce to the stress-strain relations (27) for homogeneous shells. 


CONCnJDIHG REMARKS 


In this report the principle of complementary energy is derived 
for thin plates and shells with large deflections. The plates and shells 
may he either homogeneous or of sandwich-type construction. The 
application of the principle to the derivation of the stress -displacement 
relations for homogeneous and sandwich shells as well as sandwich plates 
is given. 


Daniel Guggenheim School of Aeronautics 
College of Engineering 
New York University 

New York, N. Y. , Aiagust 15, 19^9 
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Figure 1 .- Element of cylindrical sliell with arbitrary cross section. 




